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Abstract 

We investigate the KNO scaling function of the modified negative binomial 
distribution (MNBD), because this MNBD can explain the oscillating behaviors of 
the cumulant moment observed in e + e~ annihilations and in hadronic collisions. 
By using a straightforward method and the Poisson transform we derive the KNO 
scaling function from the MNBD. The KNO form of experimental data in e + e~ 
collisions and hadronic collisions are analyzed by the KNO scaling function of the 
MNBD and that of the negative binomial distribution (NBD). The KNO scaling 
function of the MNBD describes the data as well as that of the NBD. 



1 Introduction 

Recently it has been found that cumulant moments of the multiplicity distributions both 
in e + e~ annihilations and hadronic collisions show prominent oscillatory behaviors when 
plotted as a function of their order g ||. In JI|] this behavior was attributed to the 
QCD-type of branching processes apparently taking place in those reactions. However, 
in [FJ, |3| we have shown that the same behavior of the moments emerges essentially from 
the modified negative binomial distribution (MNBD) (which actually describes the data 
much better than the negative binomial distribution (NBD) Q|). This distribution can 
be derived from the pure birth process with the initial condition given by the binomial 
distribution ||[]. 



In this paper we shall derive the KNO scaling function of the MNBD both by the 
straightforward method (i.e., proceeding to the limit of large multiplicities n and large 
average multiplicities (n) while keeping the scaling variable z = n/(n) finite and fixed) 
and by using the Poisson transform. Using this KNO scaling function we shall analyze 
the observed multiplicity distributions in e + e~ annihilations and in hadronic collisions. 

1 lt is interesting to mention here that it comes also from the concept of purely bosonic sources as 
presented recently in 
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2 KNO scaling function 



Let us remind that the MNBD is given by the following function |2|, |3], ^, |5| 

r1 4-r,i^ 

P(0) 



ri 



1 +r 2 



iV! W Vl + r 2 7 J r(j) V n / (l+r 2 > 

where iV (the number of excited hadrons) is an integer, r\ is real and r 2 > 

n = - fc a - 1 - — ^) (n) - -, 
2 V iV/ w 2' 



r 2 



I(C 2 -1 + I) W -I. (2) 



Moreover we have the generating function of the MNBD 

n(«) = X! = [1 - - 1)} N [1 - r 2 (u - 1)}~ N . (3) 

n=0 

On the other hand, the NBD has the following form 

T(n + k) ( k \ k ( k \ ~ {n+k) 
P{n) = WTWik) \W V + ¥)J 
where k > 0. Its corresponding KNO scaling function is the gamma distribution 

*{z) = ^rf kz z k ~ x - (4) 
2.1 The straightforward method 

Traditionally the KNO scaling function is derived from the multiplicity distribution P(n) 
multiplied by the corresponding mean multiplicity (ri) by going to the large multiplicity 
n and large mean multiplicity (ri) limit while keeping their ratio, z = lim n / n \-+oc n /( n ) 
fixed. In our case, starting from Eq. (0) we arrive at the following function 

^(z) = lim (n)P(n) 

n,(n)— >oo 

The parameters r[ and r 2 in Eq. (|5]) are given by 

A = i(c a -l-l)<„> 

A = \( C *- 1 + Jf)( n )> (6) 
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which are slightly different from ri, r2 given by Eq. ([!]), because (n) ^> 1. It should be 
noticed that the normalization of Eq. @ differs now from the unity, 

where the second term corresponds to the term (n)P(O) term in Eq. (jl]). 



2.2 The Poisson transform 

The KNO scaling function ^(z, t) can be also obtained by using the Poisson transform 
in which it is related to the distribution function P(n, t) by the Poisson transform 

„ , . inverse Poisson trans. T . , 

P(n,t) , *(z,t). 

Poisson trans. 

As a result we obtain in this approach that 



(n)/a J 2tt a J-oo ^ \a 



= — / e sz U(l-—t)ds, (9) 
2m Ja-ioo \n) 

where U(u,t) is the generating function of P(n,t). These equations hold also in the 
stationary function(t = 0). 

Using now the generating function Eq. (JJ), &(z) is given by the following inverse 
Laplace transform 

Then we arrive at the KNO scaling function for the MNBDQ 



2 In the actual analysis the delta function term occurring here will be approximated by the following 
expression 



8(z — e) = lim J — exp[-a(z - e) 2 /2]. 

a — >large y 27T 
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where L[ tt '(i) is the Laguerre's polynomial. In Eq. (|TTJ) the first term corresponds to 
the constant term in Eq. (10)(, or (n)P(O)). Because the normalization is now just 
the unity, in what follows we shall use this function as the KNO scaling function of the 
MNBD in the analysis of data. 



3 Analysis of experimental data 



We shall investigate now the applicability of the MNBD as presented in its KNO form 
(i.e., using Eq. ([TT|)) to the description of the observed multiplicity distributions in e + e~ 
annihilations |p|-|l3[| and in hadronic collisions [14, 15| . Table I shows obtained parameters 
C2, N and the corresponding values of the xLin • Moreover it shows other parameters r[ 
and r' 2 calculated from C2 and N according to Eq. @. 

The corresponding KNO scaling functions for our MNBD are compared with the 
KNO form of the above mentioned multiplicity distributions in Fig. la - If and Fig. 2a 
-2g. 



4 Summary and Discussion 

The KNO scaling function of the MNBD is obtained and is applied to the analysis of 
the observed multiplicity distribution in e + e~ annihilations and hadronic collisions. The 
data are also analysed by the gamma distribution which is the KNO scaling function of 
the NBD. As is seen from Table I, the Xmm vam es for the MNBD fit are almost equivalent 
to those for the NBD fit in both reactions. The result shows to be similar to the case 
of the analysis of the cumulant moments in hadronic collisions 0. On the other hand it 
should be noticed that the MNBD described the data of the cumulant moments much 
better than the NBD in e + e~ annihilations^. 



In order to know the stochastic structure of the MNBD in detail, we discuss the 
following point: The solution obtained from the branching equation of the pure birth 
process with the immigration under the initial condition of the binomial distribution 0, [5| 
is one of the extensions of both the MNBD and the NBD. Its generating function is given 

as 

(12) 



n(«) = E p < 



n)u 



[l-r x {u-l)] N [l-r 2 {u-l)] 



-k-N 



n=0 



where k is the immigration rate, and 



ri 



r-2 




-k Co - 1 - — 




k + N 



(13) 
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The MNBD is obtained by neglecting the power k in n(it). The generating function of 
the NBD is given by neglecting the power N in Eq. (0). The physical meaning of the 
immigration term k may be interpreted as a possible contribution from gluons. 
Using Eq. (pj), we have directly the KNO scaling function for Eq. fll~2| ) 




T(N + k) 




VvH^f— ^o- (u) 



This function becomes the KNO scaling function of the MNBD (Eq. (|ITD) when fc — » 0, 
and reduces to the gamma distribution (Eq. (|4])) if N = 0. 

In concrete application, we have confirmed that the discrete distribution of Eq. QT^D 
can not explain the oscillating behaviors of the cumulant moment observed in e + e~ 
annihilations and in hadronic collisions much better than the MNBD (Eq. ([[])). However, 
we are expecting at present that Eq. (|14D will become useful in analyses of data of some 
reactions at higher energies, since it has stochastic characteristics of the MNBD and the 
NBD. 
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Table caption 



Table I The parameters of the KNO scaling function of the MNBD used in the analysis 
compared with those of the NBD for the charged multiplicity in: (a) e + e~ collisions 
and (b) hadronic collisions. 



Figure captions 



Fig. 



1 The KNO form of the charged multiplicity in e + e collisions. The full circles are 
obtained from data of the following collaborations: HRS[pf, TASSOH, AMYpO 



DELPHI and OPALjTJ, U, respectively, in Fig. la- If. The broken line is the 
KNO scaling function obtained from the MNBD. 

Fig. 2 The KNO form of the charged multiplicity in hadronic collisions. The full circles 



are obtained from data of the following collaborations: ISR|14| and UA5 [|T5f | col- 
laborations, respectively, in Fig. 2a-g. The broken line is the KNO scaling function 
obtained from the MNBD. 
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Eq. dTID 


[GeV] 
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c 2 


A 


r 2 


xLn/NDF 
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29 
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1.083 ± 0.002 


0.039 ± 0.013 
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57 
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1.092 ± 0.001 
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1.092 ± 0.001 


0.012 ± 0.015 


1.957 ± 0.042 


4.5 / 25 
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Eq. §) 


[GeV] 






k 




xLn/NDF 


HRS 


29 






12.34 ± 0.30 




36.9 / 13 


TASSO 


34.8 






10.75 ± 0.12 




42.5 / 17 


AMY 


57 






11.90 ± 0.43 




19.1 / 19 


DELPHI 


91.2 






10.87 ± 0.12 




77.9 / 24 


OPAL 


91.2 






10.86 ± 0.17 




4.4 / 26 




133 






10.45 ± 0.54 




6.6 / 24 



Table I (a) 
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Eq. dlT]) 


[GeV] 


N 


c 2 


r[ 


r 2 


' v, mill / 


ISR 


30.4 


6 


1 190 ± 015 


123 ± 079 


1 880 ± 083 

-L • V7 V7 V / — 1 — V 7 • V 7 V 7 '.J 


15.6 / 15 




44.5 


6 


1.198 ± 0.005 


0.189 ± 0.030 


2.203 ± 0.038 


5.2 / 17 




52.6 


9 


1.205 ± 0.004 


0.599 ± 0.026 


2.017 ± 0.034 


4.8 / 19 




62.2 


10 


1.195 ± 0.004 


0.647 ± 0.028 


2.010 ± 0.036 


25.3 / 18 


UA5 


200 


4 


1.264 ± 0.011 


0.150 ± 0.118 


5.500 ± 0.156 


7.6 / 29 




546 


4 


1.275 ± 0.004 


0.368 ± 0.060 


7.718 ± 0.243 


54.2 / 45 
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Eq. (§ 
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Table I (b) 
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